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Abstract 

The small and negative value of the Standard Model Higgs quartic coupling at high scales can 
be understood in terms of anthropic selection on a landscape where large and negative values are 
favored: most universes have a very short-lived electroweak vacuum and typical observers are in 
universes close to the corresponding metastability boundary. We provide a simple example of 
such a landscape with a Peccei-Quinn symmetry breaking scale generated through dimensional 
transmutation and supersymmetry softly broken at an intermediate scale. Large and negative 
contributions to the Higgs quartic are typically generated on integrating out the saxion held. 
Cancellations among these contributions are forced by the anthropic requirement of a sufficiently 
long-lived electroweak vacuum, determining the multiverse distribution for the Higgs quartic in 
a similar way to that of the cosmological constant. This leads to a statistical prediction of the 
Higgs boson mass that, for a wide range of parameters, yields the observed value within the la 
statistical uncertainty of ~ 5 GeV originating from the multiverse distribution. The strong CP 
problem is solved and single-component axion dark matter is predicted, with an abundance that 
can be understood from environmental selection. A more general setting for the Higgs mass 
prediction is discussed. 
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1 Introduction 


The discovery of a perturbative Higgs at LHC, together with the absence of any signal in recent 
decades of new physics at particle accelerators or in searches for dark matter, suggests the 
possibility of a paradigm shift in Beyond Standard Model Physics: the Standard Model (SM) 
may be valid up to scales very much larger than the weak scale. The key questions then become 
very different from those associated with new physics at the TeV scale, and in this paper we 
focus on 

• The origin of an unnatural weak scale, v. 

• The origin for an unnaturally small value for CP violation in QCD, 9. 

• The origin for the small Higgs quartic coupling at very high scales, A. 

• The nature and cosmological abundance of dark matter, pm- 

The multiverse provides a framework for understanding the highly fine-tuned values of both the 
cosmological constant [1, 2], including a solution to the Why Now Problem [3, 4, 5], and the 
weak scale [6, 7]. However, it does not explain the smallness of 9. By far the most compelling 
understanding of 9 is provided by a Peccei-Quinn symmetry [8], which promotes 0 to a held. 
Furthermore the resulting axion [9, 10] can account for the observed dark matter. 

Can the multiverse provide an explanation for the A? In this paper we provide a particularly 
simple model for the axion that relies heavily on the multiverse, not just to account for the 
observed value of v and perhaps pm, but because it typically leads to a large negative value for 
the Higgs quartic coupling A at high scales. From a conventional view, where v, pm, A are to be 
understood from symmetries, this model is disastrous. But in the multiverse view, parameters 
that determine the Higgs potential must be anthropically selected to yield a universe with a 
sufficiently stable electroweak vacuum. This allows us to compute the probability distribution 
for A and obtain a successful statistical prediction for the Higgs boson mass. 

We study a supersymmetric extension of the SM that has a PQ symmetry and is valid to 
some very high energy scale M* that may be near the Planck scale. Below M* the theory 
contains no dimensionful parameters other than the scale rh of supersymmetry breaking, which 
we take to scan. The Higgs sector of the theory contains a gauge singlet held S in addition to 
two Higgs doublets with the interaction 

Ws = (1) 

where ^ is an order unity coupling. One possibility is that this theory is the MSSM with p 
replaced by S. This is not the NMSSM since S is charged under a Peccei-Quinn symmetry that 
prevents any other supersymmetric interactions of S. For weak interactions to break, the soft 
mass-squared parameter for S', Hu or Hd must scale negative, generating a new mass scale, pc, 
by dimensional transmutation. Taking pc ^ v would provide a natural understanding of the 
weak scale but, for ^ ~ 1, leads to a PQ breaking scale / ~ n. This is experimentally excluded 
so this simple supersymmetric theory has not been studied before, except for the situation that 
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^ ~ v/f <C 1 [11]. With ^ of order unity, / ~ u is likely anthropically excluded because axion 
emission prevents main sequence stars having long lifetimes [12], 

Hence, in Section 2 we study radiative PQ breaking via a dimensional transmutation induced 
by m| passing through zero at a scale fic, giving 

(S') ~ ~ /ic > V (2) 

and take 10^° GeV < fic ^ Radiative PQ breaking requires fh < /Xc, while electroweak 
symmetry breaking requires rh ^ = ^{S) ~ fic so that fh, which is taken to scan in the 

multiverse, is selected to be of order order Thus the superpartner mass scale, the PQ 
breaking scale and the dimensional transmutation scale are all comparable 

br ~ ~ /Tc- (3) 

For an alternative scheme with supersymmetry broken at an intermediate scale, and the possi¬ 
bility of an axion with / ~ m, see [13, 14], 

What is the effective theory below the scale fh / ~ /ic? The anthropic necessity of 
electroweak symmetry breaking requires a single Higgs doublet, h, to be have a mass parameter 
hne tuned to be much less than fif.- Thus this effective theory is the SM augmented by the axion 
supermultiplet. In Section 3 we study the scalar potential of this effective theory. The saxion is 
of particular interest as it has a squared mass, m^, that is one-loop suppressed compared to /r^, 
and a scalar trilinear coupling A^sh'^h, with A expected to be order fh. On integrating out the 
saxion, a negative contribution to the Higgs quartic coupling results, — A_, that is parametrically 
of order —giving a SM Higgs quartic coupling 

where A+ contains the positive D and F term contributions, and all couplings are renormalized 
at TRs- Since A is typically of order fh while rus arises only at loop level, one expects A_ S> A+, 
leading to an instability of the SM electroweak vacuum [15, 16, 17, 18, 19, 20, 21]. 

Hence, in the multiverse most universes have no stable electroweak vacuum with v -C fic- 
However, A has two contributions that are typically comparable. As the soft supersymmetry 
breaking parameters are scanned, there can be a cancellation in A yielding a stable electroweak 
vacuum with v ^ /ic, which we study in Section 4. The required cancellation in A is mild, 
about 1 in 10, very much less than the cancellation required for a light Higgs doublet. Assuming 
generic distributions for the relevant soft parameters, the distribution for A in the region of the 
cancellation takes the form 

dP{A) oc dA. (5) 

This leads to a distribution for A_ at iris that is peaked towards large values, as shown in Figure 
1. Since A+ is small and does not scan, this favors large negative values of the Higgs quartic 
A at TTis, explaining why typical observers he close to the electroweak vacuum metastability 
boundary. It is known that a multiverse distribution favoring low values of the quartic coupling 
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Figure 1: The probability distribution for A_ at scale where the total Higgs quartic coupling 
is A = A+ — A_. The shaded region is excluded as the lifetime of the electroweak vacuum is less 
than 10^*^ years, with the boundary shown for ^ = 0.3 and tan/3 = 1. 


at high scales leads to a prediction for the Higgs mass [22]; however the origin and strength of 
this distribution was unknown, leading to a large uncertainty in the prediction. Nevertheless the 
physics was clear: the multiverse distribution for the Higgs quartic at high scales explains the 
close proximity of the Higgs boson mass to the bound that follows from requiring an electroweak 
vacuum with lifetime of order 10^° years. 

In Section 4, using the distribution shown in Figure 1 following from (5), we compute a 
statistical prediction for the SM Higgs quartic coupling and therefore the Higgs boson mass. 
This prediction is shown as a function of ^ by the black curves in Figure 4, for two values each 
of Pc and tan/3. The agreement with the observed Higgs mass, shown in red, is remarkable. 
In the simple supersymmetric theory described by (1), for a wide range of pc and tan/3, the 
Higgs mass is correctly predicted within the la statistical error of about 5 GeV that originates 
from the multiverse distribution of Figure 1. 

While superpartner masses are of order m ~ pc, the detailed spectrum is highly model- 
dependent. However, the masses of the heavy Higgs doublet, the saxion and the axino are 
highly constrained, and are computed in Section 5 as a function of ^ and tan /3, in units of pc. 

In section 6 we argue that the theory yields single-component axion dark matter. If R 
parity is conserved, there is a strong upper bound on the reheat temperature after inflation to 
ensure that very large, and environmentally damaging, abundances of LSPs are not produced. 
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With rapid reheating after inflation we hnd that the PQ phase transition must occur before 
inflation, so that the dark matter abundance depends on the vacuum misalignment angle 6i. 
The observed abundance, which hxes the dimensional transmutation scale to be / 10^7^2 

GeV, can be understood from a wide range of multiverse distributions for m. For example, if 
this distribution is mild, the cost of hne-tuning the electroweak scale implies a preference for 
low values of fh and therefore /, so that the observed dark matter abundance is close to the 
minimal value that allows virialization of galactic halos [23]. 

In section 7 we show that the Higgs mass prediction of Figure 4 applies to a large class of 
theories where the SM has a completion involving a SM gauge singlet scalar at scale A > 10^° 
GeV. 

A Higgs sector described by (1) is perhaps the simplest supersymmetric axion model. One 
possibility for the PQ charges is: 1 on matter, -2 on d and 4 on S. This assignment explains 
the absence of proton decay via operators of dimension 4 and 5, is consistent with SU{5) and 
50(10) grand unihcation, and leads to R parity conservation. Neutrino masses could be Dirac, 
or the seesaw mechanism can be implemented from the vev of an addition held carrying both 
PQ and lepton charges. 

2 PQ Breaking from Dimensional Transmutation 

The model is obtained by adding a SM gauge singlet chiral superheld S to the MSSM held 
content. The interactions for S are given by the superpotential of eq. (1) together with the soft 
SUSY breaking potential 

Usoft = SHuHd + h.c. + ml |57 (6) 

At tree level the resulting lagrangian has an exact U{1) PQ global symmetry under which S and 
HuHd rotate with opposite charges, forbiding the appearance of a //-term and the corresponding 
soft SUSY breaking H/x-term. 

We dehne the various parameters and helds at renormalization scale M*, the cutoh of the 
theory, which we denote by a subscript. Both 7 and can be taken to be positive with 
no loss of generality. According to the sign of the soft parameter we can distinguish two 
cases. If < 0 then the PQ preserving vacuum is unstable. The saxion vev is stabilized 
by higher dimensional operators and (S) ~ M*. We discard this possibility. If, on the other 
hand, > 0 the PQ symmetry is unbroken at tree leveU. Once loop effects are included, 
however, the potential is modified for values of the held S M* and a PQ breaking minimum 
can develop. This occurs because the superpotential coupling in eq. (1) will generically drive 
the S soft mass negative at a scale 

fx, ~ (7) 

where gauge and top Yukawa couplings have been neglected and for illustration we take equal 
soft parameters for the scalars in Hu, H^, and S and for the trilinear interaction: = 

^An order fh value of can drive PQ symmetry breaking. Notice that in this case one expects {S) ^ 
{Hu){Hd)/fh which is unacceptable. 
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Figure 2: Left panel; phase diagram of our model for rh = 10^^ GeV, showing regions of PQ 
breaking (blue), EW symmetry breaking (red) and no symmetry breaking (gray). Right panel; 
running of the soft parameters and m| for the two benchmark points A and B. 


= y4|^. The scale /ic is a priori independent of the overall scale of supersymmetry breaking 
fh and is generated by dimensional transmutation. 

This radiative breaking of PQ symmetry is illustrated in Figure 2. The left panel shows a 
phase diagram in the (m|^,^*) plane for fh = 10^^ GeV, with other relevant parameters fixed 
at = ^u* = = hr and r/t* = 1 for the top quark Yukawa 

coupling. In the gray region neither electroweak nor PQ symmetry breaks, while in the red 
region electroweak symmetry breaks at a scale /Xc ^ fh. Hence it is the blue region, with PQ 
symmetry breaking at the dimensional transmutation scale /Xc ^ fh, that we study in this paper. 
The right panel shows running of and for the benchmark points A and B dehned in 
the right panel. 

The vev of S is computed by minimizing the full effective potential for S, which at 1-loop 
leading-log order is 

V(S) = A(,.) + m|(^)|S(,.)|" + (8) 

where is the 1-loop Goleman-Weinberg [24] potential 




2si ‘ 






log 




( 9 ) 


where the sum extends over all particles i in the model with mi{S) being their held dependent 
masses and s* their spins. Given the leading-log accuracy at which we are working, all the various 
parameters appearing in should be evaluated at a reference scale /iq such that log/r//io is 
not much larger than 1. The running parameters in eq. (8) are evaluated at the running scale 
fi according to their RG equation. The potential V is ^ independent at the order at which we 
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are working. This can be explicitly checked using 




dk 


dlog/i 


= +2rri\j + mt, = —2,^^. 

dlogp ^ 


( 10 ) 


The running cosmological constant is included to cancel the /i dependence of the held indepen¬ 
dent part of V. 

It is extremely convenient to choose the scale po to be the renormalization point pc at which 
m|(pc) = 0, so that the leading-log effective potential becomes 


1 ^( 5 ) = 
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K 


mt 


4mH log — - - + 4m;, log — - - - 8m^ log — - - 


K 


mt 


K 


(11) 


where the parameters m'jj, m\ and mp depend on S through 


m 




H,h - 2 

mp = ^|5'|. 




2 2 \ 2 
mir — mij 


+ enfisi" 


( 12 ) 

(13) 


All soft masses and couplings are now evaluated at pc so that the relevant parameter set is 

fic, mjjjfic), Mkc)- (14) 

The scale pc arises from dimensional transmutation, the soft parameters are of 

order the supersymmetry breaking scale m, and we take ^ to be order unity. Minimizing V(S) of 
(11) leads to ^(S) ~ jUc, so that the axion decay constant is given by the dimensional transmuta¬ 
tion scale / = (S) ~ /ic- This generates a supersymmetric Higgs mass parameter jj, = ^{S) ~ pc 
which is unrelated to the scale of supersymmetry breaking. If pc 3> m electroweak symmetry is 
unbroken; if pc m radiative PQ breaking fails, as the superpartners become massive before 
m| scales negative, so that either electroweak symmetry is unbroken, or PQ symmetry breaks 
at the weak scale. Hence, the viability of this model requires a special choice fh ~ pcj further¬ 
more, since f ^ supersymmetry breaking occurs many orders of magnitude above the weak 
scale, which is therefore highly fine-tuned. In the next section we take the soft supersymmetry 
breaking parameters at M* to scan in the multiverse and study the consequences of imposing an 
environmental requirement of electroweak symmetry breaking. The entire relevant parameter 
set of (14) scans; however, since ^(M*) does not scan, scans only via its pc dependence 

which is logarithmic and mild. 


3 Higgs Quartic Coupling of SM 

To obtain the SM Higgs potential we must identify the SM Higgs state and obtain the effective 
theory by integrating out super partners and heavy scalars. 



In a supersymmetric theory in which there is a hierarchy between the scale of SUSY breaking 
and the weak scale rh 3> mz, the condition for EWSB can be expressed as ~ 0, where 

Ml is the Higgs mass matrix around the EW preserving (but PQ breaking) vacuum 


M 


2 _ 
H — 


+ ml^ 

V 


+ mlj 


(15) 


where we dehned fi = Furthermore, with ^ of order unity jj, ~ /ic- As discussed above, 

the requirements of both PQ and EWSB relate the two a priori unrelated scales /ic and fh. 

On hne tuning detMl ~ 0, it is convenient to rotate from the (iP„, Hd) basis to the (h, H) 
one, where h is the massless SM Higgs doublet: 


Hu = 

sin /3/i — cos /3i7 , 

(16) 

Hd = 

cos I3h + sin I3H , 

( 17 ) 

tan /3 = 

//i2 + ml^ 

V /i2 + ml^ ■ 

(18) 


We dehned h = ia 2 h* and similarly for H. 

As shown in Figure 3, on decoupling the superpartners and expanding around the = 0, 
S = Vs vacuum, the effective theory below Hc is the SM augmented by the axion and saxion 
helds (and an axino which we discuss later). The scalar potential for this effective theory is 


777 

V (h, s) = sh'^h + \ss‘^ + \sh s'^h^h + \h (h'^h)^. 


(19) 


The mass of the saxion s is obtained by expanding the Coleman-Weinberg potential of eq. (11) 
around its minimum, leading to a value which is loop-suppressed compared to the scale of m 
and 


rnt ~ 


(47r) 




( 20 ) 


Furthermore there is a trilinear scalar interaction between s and h, A^sh'^h, where 


A = -^2 ^/i — ^ sin 2/3 


( 21 ) 


The various quartic couplings are hxed by supersymmetry and are positive sums of D- and F- 
term contributions. Matching at tree-level at the scale /ic gives 

o' 

\h = -z-cos^ 2/3 -I- ^ sin^ 2/3 (22) 

8 4 

with all couplings evaluated at the scale /ic. 

^We assume Im (S') = 0 without loss of generality. 
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Figure 3: The Effective Field Theories below M*, and the saxion mass 


As shown in Figure 3, on integrating out the saxion, the effective theory below is the SM 
augmented by the axion with scalar potential 

V{h) = \{h^hf. (23) 

Matching the theories at tree level at the boundary at rris gives 

X{ms) = Xhirus) - (24) 

2 mj 

where the last term arises from a diagram with a virtual s and two trilinear interactions and 
we ignore the running of and between /ic and The important point is that this term 
is negative-definite and very large due to the lightness of s. The natural expectation is that 
/nil ~ IbTT^, so that X{ms) <C 0 and there is no electroweak vacuum with (h) <C nig. 

However, we insist on an environmental requirement of electroweak symmetry breaking at 
a scale close to that observed. To accomplish this we scan the soft supersymmetry breaking 
parameters at M*. This allowed us to take one combination of the Higgs doublets much lighter 
than He, but this tuning alone is insufficient. If no universes have the required electroweak 
symmetry breaking then our theory is excluded. However, some universes do have the desired 
weak scale, those that have a cancellation between the two terms in (21) so that {AH/^/ml ~ 1 
(recall that fi ~ fie, which scans with the soft parameters)^. This is only a 1 in 10 fine-tune, 
much milder than that required to get a light Higgs doublet. This fine-tuning, like that for mh, 
is an inevitable consequence of our theory. The necessity of this fine-tuning has an important 
consequence: it allows us to predict the probability distribution for the Higgs quartic coupling 
and therefore the Higgs boson mass, as we discuss next. 

4 Predicting The Higgs Quartic from Vacuum Stability 

For generic values of the various soft parameters entering the Higgs sector, the SM Higgs quartic 
in eq. (24) is negative and much too large to lead to an acceptable electroweak vacuum state: 
for order one tan jS, the negative contribution coming from integrating out the saxion field is 
between one and two orders of magnitude bigger than the vacuum stability bound [15, 16, 17]. 

^ rUg also avoids having a deeper minimum of the potential (19) at (s) ^ {h) ~ A/^. 
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Let us rewrite the quartic boundary condition at of eq. (24) as 

A(ms) = A+ — A_, A_ = Aoe^. (25) 

A+ = Xh{ms) is obtained by scaling (22) from jic to Including the leading scaling from the 
top quark coupling gives 

K = ( ^ cos^^ ^ £ ln(47r/0^ • (26) 

The dimensionless quantity 

e = - = ^2 fl - ^sin2/?^ (27) 

h V 2p J 

and Ao is a large coupling of order (dvr)^, which is dehned from eq. (24). 

The soft supersymmetry breaking parameters are taken to scan at the cutoff M*. This leads 
to only a very restricted scanning of A+. There is a very mild logarithmic scanning of {g, g', yt) 
through the scale pc which scans. As tan (3 scans, it induces only mild scanning of A+ as g, g', ^ 
are roughly comparable; also, A+ has an upper bound because cos^/3 + sin^/3 = 1. 

On the other hand, A_ and e scan strongly as they depend directly on the soft parameters. 
While e is naturally of order one and can have both signs, its anthropically allowed range is 
between one and two orders of magnitude smaller, depending on so that an environmental 
cancellation is forced on the right-hand side of (27). Due to this accurate cancellation, assuming 
generic pdfs for the various soft parameters, the probability distribution for e in the anthropically 
relevant range can be accurately approximated by a flat prior, 

dP{e) oc de. (28) 

The corresponding prior distribution for A_ is 

dP(A_) oc A:^/^ dA_ . (29) 

Since the scanning of A+ is so mild, this allows us to calculate the probability distribution for 
observed values of X{ms). 

The requirement that the SM electroweak vacuum has a lifetime longer than 10^° years, 
against quantum tunneling at scale m^, leads to a constraint X{ms) > Xcr{ms), and Acr is known 
accurately [18, 19, 20, 21]. In our theory this is interpreted as an anthropic constraint on A_ 

A_ < A_|_ — Aqj" (30) 

In this region that allows observers the normalized distribution is 

^- 1/2 

~ dX.. (31) 

+ Acr 

^This same argument is used in the literature to obtain the a priori distribution for the cosmological constant 
which is then used to discuss its anthropic constraints. 


dP(A_) = - 


2 yt 
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With dP{X_) in hand, we compute the average value of the Higgs quartic and its variance 


W — 2"^+ + cr\ — 1-^+ “ (32) 

where = (A^) — (A)^. Eq. (32) shows that the average value of A gets farther and farther from 
its catastrophic value Acr and asymptote to (A) ~ |A+. 

In Fig. 4 we plot the mean and la range of X{fic) as a function of The values of eq. (32) 
are scaled up from rris to jic using SM RG equations so that all quantities in Fig. 4 refer to 
hxed values of /ic = 10^^ GeV (10^“^ GeV) in the left (right) panel. We include the 1-loop top 
contribution from (26), which is visible in the prediction at low The theoretical extrapolation 
to /ic of the low-energy Higgs quartic coupling obtained from the Higgs mass measurement [21] 
is shown by the red band. The agreement of the Higgs mass prediction (black) with data (red) 
is striking. In particular, very small values of A(/ic) in the range from —0.02 to -1-0.01 can be 
understood from values of that are order unity. The size of our Icr prediction on the Higgs 
boson mass can be estimated from the the uncertainty on the Higgs quartic running due to the 
experimental uncertainty on the measured Higgs boson mass [21]. We hnd that the la band at 
the scale fic translates into an uncertainty on the Higgs mass prediction of the order of 5 GeV. 


5 The spectrum of the model 

The three input parameters ml, and can be traded for the more physical ones det 
A and tan fd. det and A are furthermore constrained by the anthropic requirements coming 
from EWSB. To simplify the discussion of the model we will £x both of them to be zero, even 
though this is not strictly necessary for A. Having done this the Higgs sector of the model is 
completely dehned by three parameters /Xc, 'C(hc) and tan/3. In the left panel of Fig. 5 we plot 
the ratio as a function of tan/3. As anticipated in the previous section this ratio 

is 0(1) for all relevant values of tan/3. Furthermore it is independent of ^ as the parameter 
enters the lagrangian only through the combination ^S. 

In the right panel of Fig. 5 we show, as a function of tan/3, the value of the heavy Higgs, 
saxion and axino masses, normalized to the scale /Xc- The axino mass which is generated at tree 
level is tiny 

ev^ 

sin 2/3 (33) 


..tree 


m~ = 


/i 


coming from integrating out the higgsino 
contribution 

r2 


sin 2/3 


The axino mass receives a much larger one-loop 

IhP 


ml 

fj? — m\ 


In 


m: 


(34) 


This loop is analogous to the one-loop higgsino threshold contribution to the bino and wino 
masses in anomaly mediation [25, 26]. The axino typically receives a contribution of order the 
gravitino mass from higher dimensional operators [27] which we do not include here as the 
gravitino mass is model-dependent. 
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Figure 4: Statistical prediction for the SM Higgs quartic coupling X{fic) as a function of 
for /ic = 10^^ GeV (first row) and fic = 10^"^ GeV (second row), and for tan/? = 1 (left column) 
and tan/? = 2 (right column). The solid line gives the average predicted value, while the dashed 
lines show the la statistical uncertainties in the prediction. The red shaded region is allowed 
by experiment at la and 3a. 
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Figure 5: Left panel: ratio ^vs/fic as a function of tan/3. Right panel; heavy Higgs, saxion and 
axino masses normalized to the scale fic as a function of tan /3. The saxion and axino masses 
are also shown for three reference values of 

Since the RG evolution makes ^ grow in the UV, the allowed values of ^ can be constrained 
from above by the requiring it to be in the perturbative range up to the cutoff scale M*. In 
Fig. 6 we show the maximal value that ^ can attain at the scale fic imposing that < 0.3 

at the Planck scale. This maximal value depends on tan/3. In Fig. 6 we also require the top 
Yukawa coupling to be perturbative up to the Planck scale < 0.3), which eliminates 

the small tan /3 and small pc region. 

Arbitrarily small values are in principle allowed for Notice however that very small values 
of other than being difficult to explain from a theoretical point of view, induce a very slow 
running of from the cutoff scale M* down to the dimensional transmutation scale /ic- This 
implies that the UV spectrum of theories with very small ^ is rather peculiar, with m| being 
much smaller than the other superpartners soft masses. 

As the PQ symmetry of the model is broken by the vev of S the theory contains a massless 
(at tree-level) boson. We want to identify this Goldstone boson with the QGD axion. In order 
to study its interaction it is useful to parametrize the neutral components of S, Hu and Hd as 

Vs + ^ , i = u,d. (35) 

The phases (ps,i ,2 are a valid parametrization for the pseudoscalar excitations of the theory and 
Vs = {S). Only one linear superposition of the (/)’s gets a tree level mass in the theory. The 
relevant combination is easily identihed as the only dependence on the phases of the various 
held in the potential is through the A-term for ^ 

A^SH^H^ + h-c. D 2A^VuVdVs cos 
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Figure 6 : Maximal value of ^ at the scale /Xc under the assumption of a maximal allowed value 
of ^ at Mpi. In the gray region > 0.3 at the Planck scale. 

We can thus write 


G = 

cos/3(f)d - sin/3(/)„ 


(37) 

A = 

(sin/30rf + cos (3(t)u) 

cos a — sin 005 

(38) 

a = 

(sin/3(/)rf + cos/34>u] 

sin a + cos 

(39) 


where G is the linear combination eaten by the Z boson, A is the heavy pseudoscalar and a is 
the massless axion. We also have 

V V 

tan/3 = — and tan a = — sin/3 cos/3 —. (40) 

Vd Vs 

Under the following rephasing of the helds 

Hi ^ S ^ with Qi = — siii^ /3Qs, Q 2 = —cos^ l3Qs (41) 

one has 

G ^ G, A ^ A, a—)-a + fQsa (42) 

where we defined 

/ = \/2 Vs cos a — V 2 V sin a sin f5 cos /3 = \/2—— V2vs. (43) 

cos a 

In the following we £x Qs = 1. The held a is thus the canonical axion held, shifting in a 
simple way nnder a PQ transformation. We refer to Appendix D for a discussion of the leading 
interactions of the axion held in our model. 
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6 Three Scenarios for Axion Dark Matter 


6.1 Dark Matter Overview 


In general we might expect dark matter to have both axion and LSP components. In this section 
we argue that dark matter is entirely axionic, and results from three distinct multiverse scenarios 
with differing values of /. 

The superpartner mass parameters at the cutoff scale M^, scan in our theory. It is convenient 
to dehne m as the typical soft mass relevant for the hue tuning of the weak scale and, after 
marginalizing over the condition that a sufficiently stable electroweak vacuum exists, we assume 
a power-law distribution for m. The dynamically generated scale /ic ~ ^{S) also scans in the 
multiverse, and in principle is unrelated to the SUSY breaking scale. However, as extensively 
discussed in Sec. 3, the anthropic requirement of successful EWSB imposes pc ~ hr. In this 
Section we investigate the scanning of rfi in the multiverse and its consequence for dark matter. 

Although ^ does not scan it is a free parameter and, furthermore, there are other parame¬ 
ters arising from inflation that affect the dark matter abundance. An important distinction is 
whether the (last) PQ phase transition occurred before or after inflation, leading to two differ¬ 
ent cosmological axion scenarios: “Pre-Inflation” and “Post-Inflation”. The axion is not present 
during inflation if the PQ breaking scale / is lower than the Gibbons-Hawking temperature 




Hi 


E] 

V 3 Mpi 


(44) 


with Ej the energy scale of inflation. Even if this is not the case, PQ symmetry may be restored 
afterwards if / < Tmax, with Tmax the maximum temperature attained during the reheating 
process. Thus we define Pre-inflation cosmology from the requirement 


m 


/ ~ — > max (Tgh, T„ 


( 45 ) 


Throughout this section, we simplify our discussion by working in the instantaneous reheating 
limit, corresponding to a sufficiently large inflaton decay width. In such a limit there is no 
distinction between T^ax and the reheating temperature Trh, and they are obtained from the 
energy scale of inflation 

Trh = Tmax = 0.66 El . (46) 


This provides an additional parameter, which we assume does not scan. With Ei well below 
Mpi, the boundary between Pre- and Post-Inflation scenarios is at / ~ TpH- 

In the absence of R-parity violation (RPV), the lightest supersymmetric particle (LSP) is 
stable and may contribute to the energy density of the universe. In our model the LSP is likely 
to be the axino, with mass given in Eq. (33) and Eq. (34). For large SUSY breaking scales this 
may be approximated as 


mLSP ~ 


—- m sin 2/9 . 


( 47 ) 
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Figure 7: Different cosmological histories in the plane, assuming instantaneous reheat¬ 

ing after inflation and taking ^ = 0.1. Since / ~ m/^, all regions contain axion dark matter. 


For large values of the LSP mass 

"^-LSP > ^RH , (48) 

inflation dilutes any primordial LSP abundance, and the plasma produced after reheating does 
not have enough thermal energy to produce LSP particles again. In such a region the dark 
matter is only made of axions, whereas if Eq. (48) is not satished we have a two-component DM 
scenario. 

The parameter space in the (rh, Trh) plane is sketched in Fig. (7) for ^ = 0.1. In the orange 
region Eq. (45) is not satished, and therefore we are in the Post-Inhation scenario. Given 
our choice for points of the parameter space in the orange region never satisfy Eq. (48), 
thus we always have a LSP relic density. The red and blue regions are associated to the Pre- 
Inhation scenario, but only the red one has a contribution to the DM relic density from the 
LSP. Throughout the orange and red regions of Fig. (7) the LSP abundance exceeds that of 
our universe by many orders of magnitude. Furthermore this region is anthropically excluded 
because the resulting dark matter density is so large that stellar systems suffer close encounters 
[28]. Hence the only region of interest in Fig. (7) is the blue region, which has single component 
axion dark matter, and is the focus of the next sub-section.® This conclusion persists even if 

®For lower values of rh than shown in Fig. (7), the axino dark matter abundance is not catastrophic. This 
part of parameter space does not describe our universe since / ^ m/^ is too low; but it might allow observers, 
so we explore it in Section 6.3. 
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Figure 8: Contours of the axion dark matter density in the Pre-Inflation region with no LSP dark 
matter, with the initial misalignment angle hxed at 1 (0.01) in the left (right) panel. The gray 
region is excluded by the halo virialization requirement of [23], and the red region is excluded 
by extremely large LSP production during the instantaneous reheating. 


the axino mass is dominated by a contribution of order the gravitino mass, or if the gravitino is 
the LSP [27], although the size of the blue region may change. 


6.2 Pure Axion Dark Matter 


We study in more detail the blue region of Fig. (7), which has single component axion DM with 
relic density 


Pa 


po 1.67 X 



Fi{9, 


f 


1012 GeV 


1.2 


( 49 ) 


with Po the observed DM abundance. Since we are in the Pre-Inflation scenario, the initial 
misalignment angle 6i is stretched by inflation to space-time regions beyond our current horizon, 
and therefore should not be averaged. Incidentally, in the blue region we do not need to worry 
about the domain wall problem (despite the fact that our model has domain wall number 
N = 3), since inflation dilutes the associated energy density. Finally, the function Fi{6) has 
been evaluated numerically in Ref. [29]. 

In Fig. (8) we show isocontours for the axion energy density for fixed values of the initial 
misalignment angle of 0.01 and 1 in the left and right panels. The blue shaded regions have 
varying axion abundance, and we caution the reader that 6i scans as well as m. The gray region 
identifies the catastrophic virialization boundary discussed in Ref. [23], while the pink region 
has catastrophic close stellar encounters if R parity is conserved. 
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As already mentioned at the beginning of this Section, we assume a power-law prior dis¬ 
tribution for the SUSY breaking scale. Keeping in mind that we are only interested in the 
Pre-Inflation scenario we have the double differential probability distribution 


dP 

d log fh d9 


<x 0(p.- pT) 0 (p?“ - pj 


v2 ^2 _|_ p^lp 


m 


'B 


(50) 


The hrst two theta functions describe catastrophic boundaries, associated to virialization and 

close encounters, respectively. In what follows, we assume = O.SpDo and 

The factor involving v accounts for the hne-tuning to satisfy the anthropic EWSB requirement. 

The term dependent on the ratio Pa/P b is a measure factor [30, 31], and pb is the baryon energy 

density which we assume does not scan. Finally, n is an unknown parameter describing the 

distribution. 

The axion energy density pa depends on the three variables or equivalently on 

{9i, m, Trh)- Below a certain 6'min ~ Hi/i^nf) the axion density is not dependent on 9i anymore. 
The condition of no LSP relic density imposes 


m > -^Tbh 


(51) 


This gives a lower bound on the range of m where we are allowed to scan for hxed ^ and Trh- 
For example, for ^ ~ 0.1 as chosen before, we end up with the condition m > 8 x IO^Trh- For 
hxed reheating temperature, the iP/-dependent term in Eq. (49) is a decreasing function of m, 
and when the inequality in Eq. (51) is saturated its size results in 


Hi ^ ^ ^ ( Frh 

27r/ “ 27r m V 3 Mpi 1^0.66 


<0.013 ■ 

Mpi 


(52) 


Given the condition in Eq. (51), we do not want Trh to be too high, otherwise we would have 
an enormous SUSY breaking scale. To be conservative we consider Trh such that the condition 
in Eq. (51) becomes m > 10® GeV (or equivalently / > 10^® GeV). This implies that the Hj- 
dependent contribution to the axion density in Eq. (49) is completely negligible, as it becomes 
relevant for values of 9^\-a that would require a transplanckian PQ breaking scale in order to 
satisfy the virialization boundary. 

For simplicity we study two slices of the (m, 9i) scanning parameters at hxed values of 

„ _ J 10“^ hne tuned angle , . 

^ \ 1 typical value ‘ 


Fixing 9i determines the allowed range for / between the virialization and close encounters 
anthropic boundaries, and taking f = fh/^ gives 


7.9 X 10^^ GeV < - 

771 

3.6 X 10^^ GeV < j 


< 3 X 10^® GeV 

< 1.4 X 10^® GeV 


o 

1 

to 

(54) 

9i = l. 

(55) 
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Figure 9: Probability distributions for m for Oi = 1 (left panel) and 6i = 0.01 (right panel). In 
each case we show the probability distribution for three representative values of n. 

In Fig. (9) we plot the normalized probability distributions for rh for n = 2,3,4 for each of 
these slices of the multiverse and for ^ = 0.1. The examples of n = 2,3,4 shown by the red, 
green and blue curves illustrate three different multiverse scenarios for axion dark matter 

• For n = 2 the probability distribution for m is peaked towards small values, leading to 
typical observers close to the virialization boundary, as proposed in Ref. [23]. There is a 
caveat to this case, if m is further decreased by many orders of magnitude to the weak 
scale, then LSP dark matter could satisfy the virialization requirements. This is addressed 
in the next sub-section. 

• The case n = 3 features a peak in the probability distribution, as a consequence of the 
measure factor from Ref. [30]. This has the remarkable feature of explaining why the 
baryon and dark matter energy densities are comparable. 

• Finally, n = 4 has a probability distribution peaked at large values of m, so that typical 
observers are near the close encounter boundary. The proximity of our universe to this 
boundary is possible, but has not been demonstrated. 

For simplicity, above we studied hxed 9i slices of the multiverse. In fact one must study the 
probability distribution of Eq. (50) over the full range of scanning parameters (m, 6i). One again 
discovers the above three behaviors, but the corresponding values of n are affected by the larger 
scan, in particular by the possibility of 9i running to small values at large /. For a distribution 
for m that is not very steep, u < 2, one discovers that typical observers have 9i order unity and 
values of / of order 10^^ GeV, close to the virialization boundary [23]. For small values of n, it 
is the cost of electroweak symmetry breaking that prefers values of / and fh as low as allowed 
by virialization. 
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Figure 10: LSP relic density for /i > 10 TeV (left panel) and /i < 10 TeV (right panel). 


6.3 The Irrelevance of LSP Dark Matter 


We conclude this Section with the results of the freeze-out calculation for the LSP relic density, 
justifying why we did not consider this contribution in the discussion above. The LSP is an 
admixture of the neutral fermions, which are the axino and the two neutral higgsinos. For large 
values of m ~ /i, of the order of 10 TeV or larger, the LSP is mostly along the axino direction 
and has a mass as given in Eq. (47). For such large values of the superpartner mass scale, we 
can neglect the mass splitting between the two higgsino-like states and integrate out the Dirac 
held with mass /i. We hnd the effective Lagrangian 


^adHH = -sin 2/3 aa h^h + h.c. 

2/i 


(56) 


where h is the SM Higgs doublet and in this sub-section ^ This interaction mediates the 

annihilation process aa —)■ which keeps a in thermal equilibrium until freeze-out is achieved. 

The (thermally averaged) annihilation cross section reads 


rp 

{(^aa^fttWrei) = XX—^ sin^ 2/3 -- . (57) 

0 rULSP 

Since we consider large values of /i ~ ??7. we are forced to pick tan/3 ~ 1. A standard freeze-out 
calculation gives the LSP density isocontours shown in the left panel of Fig. (10). The gray area 
identihes the parameter space region where uilsp < 100 GeV and therefore the annihilation to 
SM Higgses cannot control the freeze-out. As manifest from the hgure, the LSP relic density is 
at least 6 orders of magnitude above the one we observe, therefore this region is anthropically 
excluded by close encounters. 
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This conclusion no longer holds for /i close to the weak scale. In this limit we recover the 
parameter space region of the singlet-doublet state studied in Refs. [32, 33, 34, 35], with the 
singlet mass not as a free parameter but rather one-loop suppressed. As is well known, this model 
has regions of well-tempered Xlsp where we reproduce the observed DM density. The one-loop 
suppression of the singlet mass is responsible for a very light LSP, always below 100 GeV. The 
freeze-out is controlled by annihilation to SM fermions, with (thermally averaged) cross section 




XlspXlsp" 


^.J/'^rel) — (hl3^23)" 


m 


LSP 


247rc4 


(4mLSP 


— ml 


)' + mlVl 


X] {ovf + 9%) 


OmLSP 


(58) 


The sum over the vector gvf and axial gAf couplings to the Z boson runs over all the SM 
fermions except the top quark, and it reads 


ZOb+A) =7.31 . (59) 

/ 

The suppression factor (Ri 3 V 23 )^ arises from the mass mixing (for details see Refs. [33, 34]). We 
hx ^ = 0.5, since this parameter controls the mixing and we need to mix into doublet states in 
order to have efficient annihilations. In this case tan (3 is not bounded to be close to 1 anymore, 
since all the parameters are at the weak scale. Moreover, for tan /3 = 1 the mixing factor 
( 1 / 131 / 23 )^ = 0 as a consequence of a custodial symmetry of the model in this limit, therefore 
we have to be away from this limit. We show our results in the right panel of Fig. (10), where 
we plot isocontours for the LSP relic density in the (tan f3, g) plane. We observe that there are 
anthropically viable regions, and even a line where the observed abundance is reproduced. 

Hence, regions of the multiverse with m of order the weak scale have amounts of dark matter 
from axino freeze-out that would allow observers. Furthermore, a multiverse distribution with 
n <2 would make such regions with natural electroweak symmetry breaking more probably than 
regions with much larger fh. However, there are several possibilities for avoiding this unrealistic 
situation. For example, if the reheat temperature after inflation were above about 100 TeV the 
production of LSP axinos or gravitinos by thermal scattering for freeze-in can lead to excessive 
amounts of dark matter [27]. Alternatively, for n < 2, the simple power law of m"' contained in 
Eq. (50) may apply only for a limited range of fh near the intermediate scale, breaking down 
well before TeV-scale values are reached. 


7 Generality of the Higgs Mass Prediction 

We have presented our prediction for the Higgs boson mass in a simple supersymmetric axion 
model. In fact the prediction has a much wider generality. 

Suppose the SM is valid to some high scale A > 10^° GeV, and that the theory at this scale 
involves a real SM singlet scalar s with a scalar potential 

O 

7T7 

V{h, s) = Xh (h^h)^ -|- -|- A sh^h + ■ ■ ■ (60) 
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with rris ~ A. On integrating out s, the SM Higgs quartic in the low energy effective theory at 
scale rUs becomes 


A 


A/i — 


2ml' 


( 61 ) 


An electroweak vacuum with lifetime greater than 10^° years requires A > Acr, so we must insist 
that Xh > Acr- 

Our Higgs mass prediction results from two assumptions. First, at the scale A, \h is typically 
smaller than A^/ml, so that typically A < Acr- Second, A has several contributions with at 
least one scanning in the multiverse so that there is the possibility of cancellations among these 
contributions, giving a probability distribution for A of dP oc dA at small A and the distribution 
shown in Figure 1. With these assumptions, the Higgs mass prediction is as in Figure 4, with 
fil a 1-loop factor larger than A^ and \h parameterized by ^ and tan (3 as in eq. (22). 
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A RG Equations 


In this Appendix we collect the one-loop RG equations for our model. We start from the one-loop 
running of the gauge couplings, whose evolution is governed by 


^ dfj, 


Pig) , 


Pig) 



( 62 ) 


We use the S'U(5)-normalized gi, related to the SM hypercharge by gi = ^/EJSgX We consider 
a degenerate spectrum, with all superpartners at /isusv = Ai. The solutions are 


1 _ 1 
Oiai^P) CXapT^z) 




with coefficients 


( 63 ) 


= ( 41 / 10 ,- 19 / 6 ,- 7 ), = ( 5 / 2 , 25 / 6 , 4 ). ( 64 ) 

The RG evolution of the superpotential couplings has only contributions from the wave- 
function renormalization, as a consequence of the supersymmetric non-renormalization theo- 
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[36, 37, 38]. The RG equations for the top Yukawa yt and ^ read 




^ dy 

dyt Vt 

= 


4e + 3y^-3gl--gl 


e + 6!/? - y si - 3s| - gs? 


(65) 

( 66 ) 


we 


dy 167r2 

The top Yukawa is present also when the heavy degrees of freedom are integrated out, and 
know its boundary condition at the weak scale. For the SM matter held content, the top Yukawa 
RG evolution reads [39] 


dyt 
^ dy 


SM 


Vt 

IhTT^ 


■9 . 9 9 9 17 . 

- 8^3 - 4^2 - 


The top Yukawa top is matched at the SUSY breaking scale as follows 

\/l + tan^ jd 


yt{m)\ 


MSSM 


tan f3 


yt{m)\ 


SM • 


(67) 


( 68 ) 


Among the RG equations for the soft terms, the ones for the gaugino masses are the only 
ones which can be analytically solved. The RG equations read [40] 


dMa gl ha 2M, 


dlog/i 167r2 ga 

which implies Maga'^ = const. The RG equations for the A-terms result in 


dA. 


dy IbTT^ 

dAt 1 


8^g^^ + ^A^yl + 6glM2 + -g^Mi 
^‘^^ty't + + —g^M^ + QglM2 + —glMi 


dy IbTT^ 

Finally, the RG equations for the scalar soft masses read [40] 
1 


U; = (“k+“k+Y + + 6k (™k +< + < + -4?) + 

3 

5^1 — 1 . gi 


-QglM^ - %lMl + \glS 


dm? 




d log y lOvr^ 

dm?g 1 


d log y lOvr^ 

drri^Q 1 


d log y lOvr^ 

dm?, 1 


d log y lOvr^ 


2?^ (™k + ™k + ™s + A) ~ AiA - gsPff - g9i5 
[A (™k + ™k + ™l + 7l|)] , 

2k (™k + A + A + A^) - ysIaiI - 6alMl - + gk5 

4k (mk + A + < + Al) - fglMl - 


(69) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 
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where we have introduced 


4 


2 I 2 

ttitt + m 


Q 


2mi — 2m 


d 


5 = Tr ~ m] 

B RG Analytical Solution in a Simplified Case 


(78) 


We consider a simplihed case of our model, where we neglect the top Yukawa, the top A-term, 
the gauge couplings and the soft masses. The one-loop RG equations system reads 


di 


dlog/i 
dA^ 
d log fj, 
dm% 

_ -tiu 

d log fj, 

dm\j^ 

dlogii 

dm\ 

dlogii 


2^2 


IGtt^ 

2^2 


/ 2 I 2 I 2 I /i2\ 

[mH^+mH^+ms + A^) , 


+Afj , 

{ 2 ^ 2 I 2 I a 2 \ 

+ms + AJ . 


IGtt^ 

4^2 


(79) 

(80) 
(81) 
(82) 
(83) 


This system can be completely solved analytically. We give initial conditions at the cutoff scale 
M*, and we denote couplings at the cutoff with a star subscript (e.g. ^(M*) = 

The running Yukawa coupling ^ and A-term A^ result in 


^ 4 ) = 7 




A(/r) = 


i + i.iog(4^ 

A* 


1/2 


1 + I? log ( ^ 


(84) 

(85) 


The remaining equations for the soft masses can be rewritten in a matrix form 


d 


dlog/i 






( 86 ) 


It is convenient to rotate the soft masses to another basis where the equations are decoupled. 
This is achieved by the following transformation 


( 87 ) 
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V 3 

272 

Yi 
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Yi 

2V2 
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The equations for Yi variables are decoupled 


with solutions 


dlog/i 27r2 ^ Stt^ ’ 

dY2 

-^ =0 , 

d log /i ’ 

^=0. 

d log /i ’ 




Y2{^)=C2, (! 

V'3(/i)=C^3. (! 

The integration constants Q are set by the UV boundary conditions for the soft masses 

C'l = ^ + m|, - Al) , (< 


C 3 = —^ (-mH * + 3m 

2v/2 ^ “ 


2 2 \ 


We rotate back to the original basis, and we find the solution for the running soft masses 




V61 +J,log(^) 4 j^i + |^log 




V6l + 4log(^) + 




1 + log ^ (1 + 4 log 


2 ~ ^ 




2 ^-■ 

^ 72 


( 100 ) 


C Anomalies of the PQ symmetry 

The lagrangian contains 0-terms for the various gauge group. We normalize them as 

n 03 /^a n 02 yj/a yi/a 01 D D 


( 101 ) 
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Under a chiral rotation 


( 102 ) 


where ip is one of the two component Weyl fermions of the model 

Q, m'", dP, L, e'", hi, h2, s, (103) 

the various 6*-terms transform as 

9 ^9- 2T{r^)a (104) 

where T{r^) is the Dynkin index of the representation r^. T{r^) = 1/2 for fundamental repre¬ 
sentations of SU{N) and T{r^) = if the fermion tp has charge under an abelian gauge 
group. For our model we get 

^3 6^3 — NgiplQq Y Qu<: + Qd<^)oi (105) 

02 02 ~ + Ql)c> 1 — {Ql + Q2)oi (106) 

01 01 — 2iVg ^6 (1) Qq + 3 (I) QiiC-|-3 (1) Qrfc-|-2 (1) QL + Qe’^'^O'Y (107) 

-2(2(i)'gi + 2(i)'g2)« (108) 

where Ng = 3 is the number of generations and the various color and weak-isospin multiplicities 
are easily understood. Under the assumption that the PQ symmetry is exact up to anomalies 
we obtain the following relations 

0 = Qq + Qu‘= Y Q2 = Qq + Qd‘^ + gi = Ql + Qe‘^ + gi = gi + Q2 + Qs- (loo) 

Using these eq. (105) can be simplihed 

03 03 ~ NgQso: ( 110 ) 

02 —)■ 02 — Ng^SQq Y QL)Oi Y QsOi (111) 

01 01 + NgiplQq Y QL)oi -|- (l — QsOi- (112) 

Dehning 

e2 

- 9EM^^Fg^Fg,^, = sm9wW^ Y cos 9wBg (113) 

0EM = 01 -l- 02 —t 9em Y [2 — ^Ng'j QsOi (114) 

which depends only on the combination QsC(. 

D Axion interactions 

To discuss the axion interactions at low energy it is convenient to perform a held redehnition to 
go to a basis in which all the matter helds are invariant under a PQ transformation. If is 
the charge of the held tp this is obtained by 

tp —)■ e^'^Ttp. (115) 
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Through its anomaly this rotation induces an interaction of the axion with the QCD and elec¬ 
tromagnetic dual held strengths (see Appendix C) 

+ { 1 ^, - 2) (116) 

Ag = 3 is the number of SM generations. Since a transforms in the canonical way eq. (42), 
the factor Ng identihes the domain-wall number of the model which is thus different than unity. 
Eq. (116) holds all the way down to the QCD conhnement scale. This holds because in the basis 
we are using all the matter helds are neutral under the PQ rotation and they do not induce any 
modihcation in eq. (116) as they are integrated out along the RG how. 

When the held redehnition in eq. (115) is applied to the kinetic terms of the light fermions 
the following derivative interactions are obtained 

/3dga - cos^/3 d^a ^ \ 

67^756 + — - j- (d7/.75a + S7/.75S) + -^^ 7 ^ 75 ^ 1 • (117) 

We neglected derivative couplings of the axion to the fermionic vector currents. We notice 
that no interaction of the axion with the light Higgs current — iD^WH is generated. 

Similarly to eq. (116) the couplings in eq. (117) are not renormalized as the various heavy helds 
are integrated out. 


/ sin^ (3 dgtt 
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